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We discuss radiative corrections to K e s and n e s within chiral perturbation theory with virtual 
photons and leptons. We then present the extraction of the CKM elements V U s and V u d using 
the presently available experimental input. Finally we discuss the test of CKM unitarity 
based on the present independent knowledge of Vud an d V us , and describe the prospects of 
improvement with the upcoming new B.R. measurements. 



1 Motivation 

V us and V u d are, at present, the most precisely known elements of the CKM matrix, with 
a fractional uncertainty at the level of ~ 1% and ~ 0.1% respectively. This high-precision 
information is extracted from the semileptonic transitions s —* u and d — > u, occurring in low- 
energy hadronic environments. In particular, the best determination of \V US \ is obtained from 
K — ► it iu decays {Kg$), whereas the two most stringent constraints on \V u d\ are obtained from 
super-allowed Fermi transitions (SFT), and from the neutron beta decay^. The high accuracy 
of these constraints is due to two key features: the non-renormalization of the vector current at 
zero momentum transfer in the SU(N) limit (CVC), and the Ademollo Gatto theorem^. The 
latter implies that corrections to the relevant hadronic form factors at zero momentum transfer 
are of second order in the «S?7(./V)-breaking parameters (N=2,3). 

The present level of accuracy on \ V u d\ and \V US \ is such that the contribution of \V u b\ to the 
unitarity relation 

U uu = \V ud \ 2 + \V us \ 2 + \V ub \ 2 = 1 (1) 

can be safely neglected, and the uncertainty of the first two terms is comparable. In other words, 
\V u d\ and | | lead to two independent determinations of the Cabibbo angle both at the 1% 
level. Plugging in the PDG numbers^, one finds a discrepancy between the two determinations 
at the two-sigma level. An opportunity to clarify this long-standing puzzle is offered us by the 



ongoi ng ex perimental efforts which will provide new measurements of ir e T^ and branching 
ratios I I I Precise data on 7r e 3 (B.R. ~ 10~ 8 ) will allow one to extract V u d from a theoretically 
very clean environment, while improved B.R.s will allow one to update the extraction of 
V us . In both cases, however, in order to fully take advantage of new data, an update of the 
theoretical analysis is in order. In particular, one needs to have good control over the SU (2) 
and SU(3) breaking effects, as well as on radiative corrections. Chiral Perturbation Theory 
provides a modern unified framework for studying both problems. 



2 Semileptonic decays of light mesons 

2.1 Theoretical Framework 

Weak decays of hadrons are best treated within an effective theory approach. In the case at 
hand, one starts by assuming that the Standard Model is the appropriate theory at the scale 
fi ~ 100 GeV. At scales // ~ Mw one integrates out the heavy degrees of freedom and matches 
the SM to the Fermi Theory of weak interaction, enhanced by QED and QCD, with five active 
quark flavors. The effective lagrangian governing semileptonic transitions assumes the form: 

£eff = - ^\jS cw {v) x Ojj(y) , 

Ojj = l k 7m(1 - 75) Vk X Ui 7 M (1 - 7s) Vydf . (2) 

Here k are flavor indices (running over active flavors), and V is the CKM matrix. is 
the Fermi constant governing purely leptonic charged current processes (measured in the muon 
decay). Finally, v indicates the renormalization scale. 

The effective coupling is known up to non-logarithmic terms of order aa s , and is given byH 

— 1 - "T log + 0{— 



*ew(M) = l + ^(l-Z)l°g^+°F?) (3) 



After re-summing the leading logarithms 0{{a log Mz) n , a(a s log Mz) n ] via renormalization 
group, down to the renormalization point v = M p , one finds (taking into account fermion 
thresholds): 

Sew(Afp) = 1-0230 ± 0.0003 . (4) 

The second part of the problem consists in calculating matrix elements of the type (f\Ojj\i). 
In our case, the relevant matrix element has the form: 



(ir£u\Ojj\P) ~ JiJ h 



1 + O ( a log — - , ac n 

m e 



(5) 



where Ji is the leptonic current and Jh is the hadronic current. The radiative corrections will 
in general involve large long-distance logs as well as non-logarithmic pieces (ac n ). An accurate 
calculation of semileptonic decays requires control over Jh (QCD dynamics at low energy), and 
the potentially large long-distance radiative corrections. 

The appropriate framework to address both issues is provided by yet another effective theory, 
valid in the energy range below the first QCD resonance masses: Chiral Perturbation Theory 
(ChPT)EI. Active de grees of freedom in ChPT are the eight Goldostone modes (ir, K, 77) associ- 
ated with SSB of chiral symmetry in QCD, as well as light leptons and photons. The effective 
lagrangian of ChPT is constructed according to symmetry principles, taking into account the 
SSB pattern of chiral symmetry as well as explicit symmetry breaking, due to quark masses. 
The effective action is organized by chiral power counting as an expansion in the derivatives 
(momenta) of Goldstone modes (counted as 0(p)) and the quark masses (counted as 0(p 2 ) in 



the standard approach). At a given order in p/A x (with A x ~ 4:irF n ) a finite number of local 
couplings encodes the short distance physics. Even though a full matching of QCD to ChPT 
is not available, model- independent phenomenology is pos sibl e in several cases of interest. In 
this work we need the extension of ChPT developed in Ref.^l, appropriate to deal with virtual 
photons and leptons. 



2.2 P^3 decays: matrix elements and form factors 

Neglecting for the moment radiative corrections, the invariant amplitude for the decay P(pp) — ► 
ir(p n ) l + (pi) y^iPv) factorizes in the product of leptonic and hadronic currents, the latter being 
decomposed in terms of hadronic form factors. Explicitly one has 
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(7) 



We recall that the dependence of \ Ai\ 2 on /_ is proportional to (m</Mp) 2 , thus /_ is completely 
irrelevant for the electron modes (K e 3,7r e 3). Moreover, it is customary to trade f~(t) for the 
scalar form factor fo(t), defined as: 



fo(t) = f + (t) + 



Ml - Ml 



/-(*) 



(8) 



The momentum dependence of the form factors, relevant for the integral over the phase space, 
is often described in terms of one or at most two parameters (slope and curvature at t = 0), 



f + , (t) = f + (0) (l + A +i0 ^ + A' +)0 + . . . 



0) 



In this approximation the phase space integral depends explicitly only on the slope (and curva- 
ture) parameters, and we use for it the generic notation Ip(X). With this notation, the partial 
width associated to P — > ttIv is given by: 



Tp = |K/| 2 • |/+(0)| 2 • Mp ■ I P {\) ; M P 



r 2 gU4 



(10) 



The steps necessary to extract |K/|(/ = d,s) from the measured decay rates are: (1) 
theoretical evaluation of /+(0), i.e. the SU(2) and SU(3) corrections; (2) measurement (or, if 
not available, theoretical evaluation) of the momentum dependence of f+ t o(t) (form factor slope 
and curvature parameters); (3) theoretical treatment of photonic radiative corrections. The 
short distance contribution is given by the factor <S' ew and amounts to using a different effective 
Fermi coupling, while the long-distance corrections will be discussed in the next subsection [note 
that Eq llOl is not yet general enough to account for these effects, see below]. 



2.3 Long distance radiative corrections 

A recent analysis of this issue has been performed in Ref . UJ, to which we refer for details. 
Long distance radiative corrections involve the effect of (1) loops with virtual photons and (2) 
emission of real photons: 



1. Photon loops modify the very structure of the amplitude, breaking the factorized form. 
The form factors now depend on another kinematical variable 



f±(t) - f±(t,v) 



47T 



f±(t) 



where v = (pp — Pi) 2 in P + decays and v = (p n + Pi) 2 in P° decays. The function 
T c (y, Xir) encodes universal long distance corrections, independent of strong dynamics, 
and is infrared divergent. Since the dependence on the second kinematical variable can 
be factored out (to a very good approximation), the notion of effective form factor f±(t) 
survives and proves useful in the subsequent analysis. /±(0) contains pure QCD effects 
as well as new local contributions of order a. These, together with the chiral logarithms, 
are truly structure dependent corrections, which can be described in a model independent 
way within the ChPT approach. 

2. One has to consider how radiation of real photons affects the various observables (e.g. 
Dalitz Plot density, spectra, branching ratios). For the purpose of extracting we 
need to assess the effect of real photon emission on the partial widths. As is well known, 
a given experiment measures an inclusive sum of the parent mode and radiative modes: 

dr obs = dr{p e3 ) + dT(p £3y ) + ... 

From the theoretical point of view, only such an inclusive sum is free of infrared sin- 
gularities. At the precision we aim to work at, a meaningful comparison of theory and 
experiment can be done only once a clear definition of the inclusive observable is given^l 

In summary, all the universal QED effects, due to both virtual photons [T c (v, Xir)} and real 
photons [dT(K£3y)], can be combined to produce a correction to the phase space factor in the 
expression for the decay width: 



Jp(A) — I P (A,a) = I P (X,0) 



1 + A/ P (A) 



This term comes, in principle, with no theoretical uncertainty. The structure dependent elec- 
tromagnetic corrections, as well as SU(3) breaking corrections are in the form factor f+(t), 
where all of the theoretical uncertainty concentrates. At zero momentum transfer, this can be 
conveniently written as 

U(0) = f + (0) (1 + S SU (2) + Sem) 
Explicit expressions for f±(t) at NLO in ChPT, including the SU(2) and EM corrections, can 



be found in Refs.H^ 



3 V ud from 7r e3 decay 

Theoretically this process is very attractive because it shares the advantages of both Fermi 
transitions (pure vector transition, no axial-vector contribution) and neutron /3-decay (no nu- 
clear structure dependent radiative corrections). The difficulty here lies in the extremely small 
branching ratio, of order 10~ 8 . With the notations introduced in the previous sections, V U( i is 
related to the partial width by the following relation: 
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(11) 



Table 1: Summary of isospin-breaking factors. 





$SU(2) (%) 


«5em(%) 


AJ(A)(%) 


^ e3 

K% 
K% 


2.4 ± 0.2 


2.4 ± 0.2 



0.32 ± 0.16 
0.46 ± 0.08 
0.006 ± 0.16 
0.15 ± 0.08 


-1.27 

-0.32 
± 1.0 * 
1.7 ± 1.0 * 



Explicit calculation of the radiative corrections and SU(2) breaking leads tcM^l 

$SU(2) ~ !0~ 5 , %2 P 2 = (0.46 ± 0.05)% , AJ 7r (A) = 0.1% , (12) 

with a total effect of radiative corrections consistent with previous estimates'^! 
The present experimental precision for the branching ratio of the pionic beta decay cannot 
compete yet with the very small theoretical uncertainties of SFT and neutron beta decay: using 
the latest PDG valued BR = (1.025 ± 0.034) x 10~ 8 , we find 

\V ud \ = 0.9675 ± 0.0160 exp ± 0.0005 th = 0.9675 ± 0.0161. (13) 

However, a substantial improvement of the experimental accuracy is to be expected in the 
near future. Inserting the present preliminary result obtained by the PIBETA Collaboration^, 
BR = (1.044 ± 0.007(stat.) ± 0.009(syst.)) x 10~ 8 , we find 

\V ud \ = 0.9765 ± 0.0056 exp ± 0.0005 th = 0.9765 ± 0.0056 , (14) 

where the error should be reduced by about a factor of 3 at the end of the experiment. 



4 V us from K e ^ decays 



Each of the four widths is related to V us by a relation similar to Ea llll modulo an obvious 
change needed to take into account SU(3) breaking effects. The K(z widths thus allow us 
to obtain four determinations f |V^ S | • f+° n (0) which are independent, up to the small 
correlations of theoretical uncertainties from isospin-breaking corrections <5em an d 5su(2) (almost 
negligible at present, see Table The master formula for a combined analysis of these modes 
is: 



\Vus\-f. 
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(15) 



where the index n runs over the four modes (n = K^, K^, K%, K^) . In writing the above 

equation, we have chosen to "normalize" the form factors for the various modes to /£ (0), 
evaluated in absence of electromagnetic corrections. Differences between the various modes are 
due to isospin breaking effects, both of strong {5 l SU f 2 \ ) and electromagnetic (#em) origin, which 

have been evaluated at 0(ep 2 ) in the chiral expansion^! 



ft(0) 



(0) 



1 + ^SU(2) + ^Em) 



(16) 



Let us now discuss the theoretical input needed in Ea ll5l In the Tabled we report estimates 
of the isospin-breaking parameters^] a The uncertainty on 5^ M and S^j^^ is due to incomplete 



a These are based on Ref.El. 
within a different framework. 



SU{2) 

Another calculation^^ for the K^ 3 mode leads to consistent results, though 



Table 2: branching ratios (BR) and slopes from PDG fits. The lifetimes used as input are: t k ± — 

(1.2384 ± 0.0024) x 10~ 8 s and t Kl = (5.17 ± 0.04) x 10 -8 s. 



Mode 


BR (%) 


A+ 


Ao 


K% 
K% 


4.87 ± 0.06 
38.79 ± 0.27 
3.27 ± 0.06 
27.18 ± 0.25 


0.0278 ± 0.0019 
0.0291 ± 0.0018 
0.033 ± 0.010 
0.033 ± 0.005 


0.004 ± 0.009 
0.027 ± 0.006 



knowledge of the relevant NLO ChPT couplings. This is not the dominant uncertainty in the 
final analysis. T he pha se space corrections refer to the definition of photon-inclusive width 
given by Ginsberg ^ Although our calculation of phase space corrections for muonic modes 
is in progress, in order to include these modes in the phenomenological analysis, we can use 
the estimates AI n ,n = K~^ 3 ,K® 3 obtained by Ginsberg Due to potentially large model 
uncertainty (e.g. dependence on UV cutoff) we assign an error bar of ±1% to these entries. 

The expansion of fF n (0) in the quark masses has been analyzed up to the next-to-next- 
to-leading order 1 ^. At this level of accuracy we write 

/f--(0) = l + / 2 + / 4 + 0(p 6 ), (17) 

where the identity / = 1 follows from current conservation in the chiral limit. Because of the 
Ademollo-Gatto theorem^, local terms are not allowed to contribute to /2, and this implies a 
determination which is practically free of uncertainties: 

h = -0.023 . (18) 

As for the situation is much less clear, because low energy constants (LECs) of the p 6 La- 
grangian can now contribute. Various estimates of the size of f& have been given, and for a 
mini -rev iew we refer to Ref.^. In our numerical analysis we use the result of Leutwyler and 
Roos^^, based on a model-independent parameterization of the asymmetry between kaon and 
pion wave functions: 

U = -0.016 ± 0.008 . (19) 

We take this result with the understanding that, e.g., a value of f± two sigmas away from the 
central value in Ea ll9l is not strictly forbidden, but rather unlikely. 

From the point of view of the pure chiral expansion, the only parameter-free prediction which 
one can make for con cerns the chiral logs. A first step in this direction was made by Bijnens, 
Ecker and ColangelJ^, who calculated the double chiral logs contribution to this quantity. The 
size of this term, however, depends on the renormalization scale (i. By varying the latter within 
a reasonable range, the numeri cal estimate | /4 1 chiral logs < 0.5% was obtained. 

Recently two groups have performed a full calculation of f± to 0(p e ), which besides 
the double chiral logs contains single ones and polynomial contributions. The latter contain 
LECs of the p 6 Lagrangian, whose v alue is basically unknown, and make a numerical estimate 
difficult. However, a key observation^! is that the relevant local couplings turn out to govern 
the slope and curvature of the scalar form factor /q ""(i) at t = 0. Therefore, at least in principle, 
such couplings can be related to observables, through a measurement of the the t dependence 
of /rf' 7r (i) in K ^3 decays or through dispersive analyses of the scalar form factor. Finally, let us 
note that lattice calculations of /+ (0) are certainly called for in order to improve our confidence 
in the range used for 

Using the experimental input reported in Table [21 an average over the electronic modes (the 
muonic modes are irrelevant in the average, given the present uncertainties) leads to: 

\V US \ ■ /+°M0) =0.2115 ±0.0015 . (20) 
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Figure 1: \V US \ ■ /+ T (0) from the four Kiz modes (and average over the electronic modes). 



The error has been multiplied by a scale factor S = 1.54, as defined by the PDG^, in order to 
take into account the slight inconsistency between K^ 3 and K® 3 results. Finally, if we use the 
Leutwyler-Roos estimate of / 4 , or /* w ~ (0) = 0.961 ± 0.008, we obtain 

\Vus 



0.2201 ± 0.0016 exp ± 0.0018th(/ 4 ) 



0.2201 ± 0.0024 . 



(21) 



5 V u d, V us , and CKM unitarity 

Given the tiny size of V u b, the CKM unitarity constraint of Eq^ provides a relation between 
V u d and V us . Errors on the direct determinations of Kd and V us are at the level of 0.1% and 
1%. Explicitly - one has: 

I Kid | direct = 0.9739 ± 0.0005 (Dominated by n-j3 and nuclear SFT) ; 
iKsldirect = 0.2201 ±0.0024 (Dominated by K e3 ) . (22) 

CKM unitarity implies an independent determination of V us (from V uc [) at the 1% level, 

IKwUit = 0.2269 ± 0.0021 (23) 

to be compared with the direct determination in Eg 1221 The ~ 2a discrepancy between these two 
determinations could be attributed to an underestimate of i) theoretical uncertainties involved in 
I Kid | direct ] h) theoretical uncertainty in f± ; hi) experimental errors in B.R.s. (as probably 
hinted by the marginal consistency of / + (0)|Kxs| from and K^). At the moment, in absence 
of a clear indication of which of the errors is underestimated, it's probably best to treat the 
two determinations on the same footing, and introduce the PDG scale factor in the final error. 
Following this procedure one finds 

iKslunit.+A^ = °- 2240 ± °- 0034 • ( 24 ) 



6 Final remarks 



The pion beta decay (^3) provides in principle a unique test on the existing extractions of Kd 
from nuclear SFT and neutron beta decay. Theory being extremely clean, one has to wait for 
the final results from PIBETA - to address the impact of this channel on the determination of 

Kxd- 



As for Kiz decays, our theoretical understanding of radiative corrections is under control 
(implying an error on V us of about 0.2 % ), while improvement is needed in the calculation 
of SU (3) breaking effects. The experimental situation will soon improve thanks to the high- 
statistics measurements of K/% widths expected from recently-completed or ongoing experiments, 
such as BNL-E865 5 , KLOeIS and NA48^. It should be noted how a single present-day mea- 
surement allows us to extract |V^ S | at the ~ 1% level [same as in Eq|^. For example, using the 
final result from BNL-E865EI BR(J£+r ,) = (5.13 ± 0.02 stat . ± 0.09 syst . ± 0.04 norm .)%, and the 
appropriate phase space correction (AI K + = (—1.27 + 0.5)%), one obtains 

e3 

\V us \e865 = 0.2272 ± 0.0030 . (25) 

Although this result is in perfect agreement with what is expected from unitarity, we stress 
that when combining BNL-E865 with other existing results one finds: 1) internal inconsistencies 
within ifij measurements; 2) the discrepancy with K® 3 is worsened; 3) in a combined analysis, 
the problem with unitarity persists at the ~ 2a level. In our opinion, a meaningful improvement 
will be possible once the database is fully updated, including not only BNL-E865 but also 
the results from KLOE and NA48. 
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